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ABSTEACT 


The  pilot  wind  tunnel  is  a  variable  geometry  wall  transonic 
tunnel.  It  is  constructed  to  use  air  from  an  existing  bottle 
farm  and  blow  down  to  atmosphere.  The  report  presents  the 
criteria  for  selection  of  control  components  and  discusses 
in  detail  the  design  of  the  airflow  control  systems. 
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FOREWARD 


The  work  described  herein  is  part  of  an  on-going  project  to 
verify  a  concept  for  a  wind  tunnel  flexible  wall  test  section 
in  a  pilot  scale  facility.  This  effort,  including  the  actual 
pilot  tunnel  which  is  approaching  operational  status,  is  being 
carried  out  at  the  Air  Force  Flight  Dynamics  Laboratory,  Wright- 
Patterson  Air  Force  Base,  Ohio  45433.  The  work  is  an  element  to 
in-house  Work  Unit  Number  24041303  titled  "Develop  New  Technology 
Concepts  for  Aerodynamic  Ground  Testing"  under  Task  Number  240413 
titled  "Aerodynamic  Ground  Test  Technology”. 

Aerodynamic  guidance  and  criteria  for  this  effort  are  being 
supplied  by  Dr,  Donald  J.  Harney,  Assistant  for  Simulation  to  the 

t  '  * 

Aeromechanics  Division  Chief* 

The  valuable  assistance  of  many  technicians-  and  several  engi¬ 
neers,  especially  Mr*  Alan  G.  Blore,  in  assembling  this  tunnel  is 
greatly  appreciated. 

The  author  wishes  to  thank  Ms.  Elona  Beans  for  typing  the 
report  and  several  illustrators,  especially  SSgt  Nicholas  J. 
Novogradac,  for  preparing  the  artwork. 
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I.  INTRODUCTION 


In  late  1973  and  early  1974  plans  were  being  formulated  in  the  Air  Force 
Flight  Dynamics  Laboratory  (AFFDL)  for  constructing  a  wind  tunnel  with  a  fifty 
inch  square  test  section  with  variable  geometry  walls.  The  variable  geometry 
of  the  walls  would  provide  an  adaptive  test  section  which  would  allow  testing 
larger  models  in  a  smaller  tunnel  than  possible  in  a  comparable  sized  conven¬ 
tional  tunnel.  The  facility  would  provide  support  for  in-house  test  programs 
when  time  in  other,  larger  tunnels  was  not  available,  and  provide  this  support 
at  a  lower  cost.  These  economies  would  also  permit  a  more  thorough  job  of 

testing,  which  is  important  in  a  flight  regime  where  discontinuities  exists  in 
modeling  theories. 

In  order  to  assure  that  the  design  of  the  variable  geometry  is  accurate. 
Planning  was  initiated  for  in-house  construction  of  a  blowdown  pilot  facility 
with  a  nine  inch  square  test  section.  The  pilot  facility  would  provide  a  means 
of  testing  the  variable  geometry  concept.  It  would  also  provide  a  means  of 
determining  the  degree  to  which  the  hardware  realization  would  have  to  match 
theoretical  predictions  of  wall  curvature  requirements.  Operating  experi- 
with  the  pilot  facility  should  permit  relaxing  the  curvature  requirements 
and  result  in  a  test  section  of  less  complexity. 

This  report  does  not  cover  the  aerodynamic  design  of  the  test  section  wall 

ture.  (Harney,  1973).  The  report  is  intended  to  document  the  design  of 
the  test  leg  air  controls. 
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Design  of  test  section  is  documented  in  an  AFFDL  Technical  Report  (TR) 
preparation.  ' 


II.  DESCRIPTION  OF  FACILITY 


A  brief  description  of  the  tunnel  airflow  circuit  and  the  test  section  is 
included  herein.  Test  section  details  are  described  in  this  report.  Airflow 
circuit  details  are  described  in  a  Technical  Report  in  preparation.  A  schematic 
diagram  is  shown  in  Figure  1. 

A.  TUNNEL  CIRCUIT 

Air  Storage.  The  tunnel  will  operate  from  the  same  high  pressure  bottle 
farm  as  the  AFFDL  Mach  3  High  Reynolds  Number  Facility  (Fiore,  1975).  This 
15,000  Ft  of  3,000  psi  air  storage  is  shared  by  a  number  of  wind  tunnels  and 
test  rigs  in  the  Bldg  450  complex.  Due  to  line  losses  and  other  air  usage,  the 
useable  air  for  this  facility  is  approximately  15,000  pounds.  The  recharge  rate 
varies  from  15  to  2  lbm/ sec  depending  on  how  many  of  four  available  reciprocating 
compressors  are  on  the  line.  The  air  consumption/recovery  ratio  is  approximately 
100  to  1.  A  two  minute  run  would  require  more  than  three  hours  of  pumping  time 
during  which  some  facilities  could  not  run.  Therefore  the  wind  tunnel  design 
must  maximize  productivity  for  very  short  runs.  The  air  circuit  controls  are 
designed  to  accommodate  runs  of  twenty  seconds  or  less. 

2 •  Supply  Piping.  Supply  piping  to  the  test  leg  is  shown  in  a  schematic  of 
the  facility  Figure  1.  The  perforated  plate  will  act  as  a  flow  distributive 
throttling  device  to  remove  the  control  valve  signature  which  sometimes  affects 
wind  tunnel  data  (Lowe  and  Cumming,  1970) . 

The  piping  downstream  of  the  regulating  valves  is  protected  from  over 
pressure  by  a  rupture  disk. 
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FIGURE  1  -  Schema!: ic  of  Facility 


3.  Test  Leg.  The  test  leg  consists  of  a  settling  chamber,  a  three  dimensional 
butterfly  contraction  section,  a  two  dimensional  nozzle,  a  test  section,  a  rectan¬ 
gular  to  round  transition  section  and  a  hydraulically  operated  butterfly  style  con¬ 
trol  valve.  This  control  valve  provides  a  means  of  choking  the  flow  downstream  of 
the  nozzle  for  setting  subsonic  Mach  numbers. 

The  settling  chamber  has  a  cross  sectional  area  of  thirty-one  times  that  of 
the  test  section.  It  contains  a  conical  baffle  and  screens  to  reduce  turbulence  in 
the  flow. 

Since  numerous  runs  are  planned  with  P0  close  to  the  design  pressure  limit, 
two  safety  devices  are  incorporated  in  the  settling  chamber.  The  low  pressure 
safety  device  is  a  relief  valve  which  vents  the  settling  chamber  to  exhaust  line 
and  silencer.  Should  the  pressure  continue  to  rise;  a  rupture  disk,  capable  of 
handling  a  much  larger  mass  flow,  vents  the  settling  chamber  to  atmosphere. 

4.  Exhaust  Piping.  The  butterfly  control  valve  is  connected  to  a  silencer 
with  two  foot  diameter  low  pressure  piping.  For  personnel  safety,  the  stream 
exhausts  vertically. 

The  plenum  areas  on  each  side  of  the  test  section  are  connected  to  a  single 
manually  operated  butterfly  style  control  valve.  The  control  valve  is  used  to 
modulate  the  plenum  flow.  The  plenum  flow  can  be  returned  to  the  stream  at  the 
silencer  inlet  when  sufficient  differential  pressure  exists.  Lack  of  sufficient 
pressure  drop  would  necessitate  connecting  the  plenum  line  to  existing  vacuum 
source  in  the  testing  complex. 
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B*  TEST  SECTION.  The  test  section  wall  deflection  capability  is  shown  diagram- 
matically  in  Figure  2  and  is  discussed  further  in  Section  III.  The  figure  scale 
is  not  that  of  the  pilot  tunnel,  and  the  deflection  scale  is  exaggerated.  The 
walls  will  be  contoured  by  jacking  to  provide  minimum  wall  interference  in  the 
model  area  at  transonic  speeds  (Harney,  1973). 

* 

The  upper  and  lower  walls  of  the  test  section  are  each  composed  of  nine  very 
flexible  rods,  or  elements.  The  elements  begin  at  the  entrance  of  the  two  dimen¬ 
sional  nozzle  and  are  carried  through  the  test  section  to  the  transition  section. 
A  cross  schematic  of  the  test  section  is  shown  in  Figure  3.  The  back  side  of  the 
elements,  away  from  the  tunnel  centerline,  are  sealed  with  triangular  shaped 
strips  which  can  be  moved  into  close  proximity  with  the  elements.  The  seal  strips 
are  spring  loaded  against  the  elements  to  act  as  relief  devices  to  protect  against 
a  high  differential  pressure  on  tunnel  start.  The  seals  can  also  be  moved  away 
from  the  elements  to  create  a  flow  path  into  the  plenum  space  behind  the  elements. 

A  crescent  style  model  support  is  located  at  the  downstream  end  of  the  test 
section.  The  upper  and  lower  walls  diverge  at  this  point  to  compensate  for  the 
blockage  of  the  model  support. 

At  the  model  support  area  the  seals  pivot  about  a  hinge  at  the  downstream  end 
of  the  elements  and  act  as  partial  flaps.  This  provides  a  means  for  partially 
equalizing  pressure  across  the  elements.  This  also  can  provide  a  means  of  re¬ 
turning  plenum  flow  to  the  stream. 

Each  element  has  ten  Individual  jack  points  in  the  model  area  which  are 
electro-mechanically  operated.  The  walls  will  be  deflected  by  jacks  to  provide 
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FIGURE  2  -  Example  Test  Section  Wall  Contouring 
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FIGURE  3  -  Schematic  Cross  Section  Through  Test  Section 
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both  axial  and  lateral  contouring.  Between  the  nozzle  throat  and  the  model  area 
the  elements  are  connected  together  with  three  crossbeams  per  wall  at  manual 
jack  stations.  The  manual  jacks  will  be  used  to  form  a  supersonic  wall  contour 
to  develop  Mach  numbers  up  to  1.25. 
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III.  TUNNEL  CIRCUIT  DETAILS 


A.  AIR  FLOW  COMPONENTS  AND  CALCULATIONS 


1.  c-nononts  and  Abroach.  Computations  o£  flow  and  pressure  drops  for  the 
tunnel  begin  at  the  settling  chamber  and  npzrle  and  proceed  upstream  and  downstream 
from  that  point.  The  maximum  design  settling  chamber  pressure,  fixed  by  mecha 
considerations  at  sixty  psia,  and  the  nozzle  area,  fixed  by  the  aerodynamicist  at 
nominally  ei8hty-one  square  Inches,  determine  the  maximum  tunnel  flow.  The  maxi- 
mtm  tunnel  flow  is  then  available  for  the  sizing  P‘P*  end  valves.  The  tempera¬ 
ture  Is  not  controlled  and  is  assumed  to  be  standard,  520  K. 

Proceeding  in  the  upstream  direction,  the  inlet  pressure  to  the  perforated 
plate  is  computed  first.  This  pressure  becomes  the  stagnation  pressure  control 

valve  exhaust  pressure. 


The  size  of  the  control  valve  and  of  components  upstream  of  the  valve  were 
determined  based  on  this  pressure  and  the  supply  pressure.  These  components  were 
originally  sized  for  an  installation  of  the  tunnel  in  a  different  building. 


Therefore  components  and  pipe  sizes  were  selected  to  take  maximum  advantage  of 
piping  in  place.  The  design  also  considered  the  need  for  a  large  flow,  regulated 


air  pressure  source  for  another  facility.  Based  on  these  constraints  two  regulat¬ 
ing  valves  operating  in  parallel  were  selected.  These  valves  reduce  the  supply 
pressure  from  3000  psi  to  1100  psi  or  less.  *  large  size  control  valve  which  would 
pass  the  required  flow  with  little  pressure  drop  was  selected.  A  low  pressure  air 
actuator  was  selected  for  this  balanced  control  valve  with  the  intention  of  pre¬ 
setting  the  valve  and  only  fine  tuning  the  control  with  flow  on.  However,  since 
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other  facilities  use  the  stored  air  in  the  new  location  which  was  chosen,  this 
valve  is  now  being  used  to  bring  the  tunnel  on  line  directly  and  rapidly. 

An  economical  design  for  the  present  installation,  starting  from  scratch, 
would  have  dictated  the  use  of  a  high  pressure  valve  to  control  the  flow  with 
only  one  stage  of  control.  However  the  use  of  suitable  components  which  are 
already  purchased  is  necessary  for  economy,  and  the  original  vhlves  were  retained. 
Since  the  piping  from  the  regulators  to  the  control  valve  is  much  shorter,  the  line 
losses  are  much  lower  and  pressure  fluctuations  for  flow  versus  no  flow  are  much 
less,  control  with  the  valve  is  easier  than  it  was  for  the  original  installation. 
The  valve  actuator  hhs  also  been  replaced  with  a  hydraulic  cylinder  and  servo¬ 
valve  to  provide  a  quick  start  capability  which  will  conserve  air. 

With  the  control  valve  size  fixed,  its  inlet  pressure  is  determined  by 
the  desired  valve  opening  at  design  conditions.  A  value  of  one-half  open  was 
selected  to  permit  an  adequate  reserve  of  extra  flotf  capacity  for  response. 

'  1118  Pressure  setting  for  the  regulating  valves  is  equal  to  the  sum  of 

the  control  valve  inlet  pressure  and  the  line  losses  between  the  regulators  and 
the  control  valve. 

The  line  connecting  the  regulators  to  the  tank  farm  is  existing  and  ques¬ 
tionable  in  size.  Its  adequacy  will  be  based  on  flow  rate,  i.e.,  the  time  re¬ 
quired  to  reduce  the  source  pressure  below  useable  values.  Should  it  prove  in¬ 
adequate  a  second  line;  in  parallel,  wLll  be  added. 

Proceeding  downstream  from  the  settling  chamber;  at  sixty  psia  total 

;  . 

pressure,  the  flow  will  exhaust  to  atmosphere.  With  a  pressure  ratio  of  four, 
sonic  flow  will  exist  at  some  point  in  this  section.  To  provide  subsonic  flow 


at  sixty  psia  total  pressure  in  the  test  section,  the  sonic  flow  section  oust  be 
located  downstream  of  the  test  section.  A  butterfly  control  valve  is  used  to  choke 
the  flow  downstream  of  the  test  section.  This  style  valve  was  selected  because 
coarse  valve  configurations,  such  as  gate  valves,  can  show  a  dependence  of  sub* 
sonic  Mach  numbers  on  stagnation  pressure,  i.e.,  the  valve  flow  coefficient  varies 
somewhat  with  total  pressure  (Grunnet,  1974).  Additionally  the  butterfly  valve 
must  not  be  choked  with  Mach  1.25  flow  in  the  test  section.  This  means  that 
when  the  valve  is  wide  open  the  flow  area  must  be  greater  than  the  eighty-one 
square  inches  nozzle  area  or  the  expanded  Mach  1.25  flow  area  of  84.8  square  inches. 

The  valve  was  selected  by  examining  the  properties  of  butterfly  valves 
ranging  in  size  from  twelve  to  sixteen  inches.  The  flow  through  the  valve  in 
the  design  choke  condition  creates  large  forces  on  the  disk.  Holding  the  valve 
open  against  these  forces  requires  a  heavy  shaft  to  provide  the  torque.  The  heavy 
shaft  means  a  thick  disk  with  less  flow  area  in  the  full  open  position.  It  was 
necessary  to  use  a  sixteen  inch  valve  with  a  class  five  shaft  (two  inch  diameter) . 

The  full  open  area  is  approximately  84.7  square  inches. 

A  silencer  was  incorporated  to  reduce  the  exhaust  noise  to  below  safe 
sound  levels.  The  silencer  is  rated  at  an  attenuation  of  48  DB  for  noise  at  134 
DB  and  100  Hz.  Effectiveness  diminishes  with  frequency  to  5  DB  attenuation  at  31.5 
Hz.  The  silencer  is  rated  at  a  pressure  drop  of  0.3  psi  at  102  lb/sec  flow  rate. 

The  two  foot  diameter  connecting  line  also  has  a  low  pressure  drop. 

The  auxiliary  flow  lines  must  also  be  adequately  sized.  Most  of  these  will  be 
discussed  under  safety  relief  flow.  Plenum  flow  results  from  moving  the  seal  strips 
away  from  the  elements  creating  a  flow  passage  through  which  a  low  pressure  down¬ 
stream  will  cause  flow.  The  maximum  flow  desired  is  only  12.5%  of  the  main  line 
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flow,  and  an  existing  ten  inch  butterfly  valve  was  found  to  be  adequate  for  throt 
tling  the  flow.  Manual  presetting  of  the  valve  is  presently  planned  but  future 
operations  may  be  under  servo  control. 

A  jury  rig  set  up  was  fabricated  and  tested  to  verify  the  airflow  quanti¬ 
ties.  The  flow  versus  the  seal  notion,  normal  to  the  test  section  rod  wall,  is 
initially  nonlinear  as  the  seal  begins  to  move  away  from  the  rods.  After  about 
10%  travel  the  characteristics  become  fairly  linear  until  about  95%  of  maximum 
flow  is  reached.  At  maximum  flow  the  spacing  between  the  rods  is  the  limitirig 
factor  rather  than  the  seal  travel. 

'  '  1  ■  i"  I  !  :  ■  .*  1.  •  t  ,  ♦  , 

2’  .E9uatlon8  and  Calculations.  Flow  through  pipes  is  assumed  to  follow  the 
Fanno  Line  (VanWylen,  1966)  for  compressible  adiabatic  flow  with  friction.  This 
method  of  analysis  was  chosen  for  pipe  sizing  because  it  considers  that  compres¬ 
sible  flow  in  a  duct  can  enter  at  subsonic  velocity  and  choke  before  the  end  of 
the  duct  is  reached. 

*  '  #  « 

Flow  through  valves  which  have  large  pressure  drop  and  low  recovery  was 
assumed  to  follow  the  standard  valve  equations  which  employ  a  capacity  factor,  Cv 

The  Fisher  Controls  equation  for  high  recovery  valves  was  used  for  the  butterfly 
style  subsonic  choke  valve. 

Flow  through  the  nozzle  is  assumed  to  result  from  still  air  in  the  set¬ 
tling  chamber  expanding  isentropically  through  the  nozzle.  This  is  a  good  assump¬ 
tion  with  a  well  made  nozzle  and  with  the  area  ratio  of  settling  chamber  to  test 
section  of  thirty  to  one. 
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(a)  Supply  line  flow  -  The  regulator  inlet  pressure  which  results  from  the 
four  inch  supply  line  losses  is  shown  in  Figure  4.  The  line  is  approximately  300 
feet  long,  has  an  inner  diameter  of  3.5  inches,  with  six  ninety  degree  elbows  and 
two  gate  valves.  The  wide  open  regulator  flow,  no  regulation,  is  shown  in  Figure  5. 
Data  for  the  latter  figure  was  calculated  using  characteristics  provided  by  the 
valve  manufacturer  (Grove,  1969)  and  by  extrapolating  for  subsonic  flow. 

Figure  5  shows  that  the  minimum  regulator  inlet  pressure  to  obtain 
the  maximum  continuous  tunnel  flow  of  115  lb/sec  is  1300  psi.  Figure  4  shows  that 
the  minimum  tank  farm  pressure  to  achieve  1300  psi  regulator  inlet  pressure  at 
that  mass  flow  is  2100  psi.  With  a  beginning  pressure  of  2600  psi,  available 
storage  of  8250  cubic  feet,  the  run  time  at  115  lb/sec  is  134  seconds  assuming 
adiabatic  expansion  of  remaining  air.  This  is  adequate  since  the  desired  run  time 
is  ten  seconds  to  sixty  seconds. 

The  pressure  drop  in  the  two  four  inch  regulated  pressure  lines  is  20 
psi  at  115  lb/sec  and  700  psi  inlet  pressure,  the  regulator  setting. 

(b)  PQ  controlling  flow  -  The  equations  for  sonic  and  subsonic  flow  through 
the  control  valve  are  (Beard,  1969): 

sonic 

gL  _2L_  p 
v^GT  2.5  1 


m  =  (.00002122) 


FLOW,  #/SEC  REGULATOR  INLET 


FIGURE  4  -  High  Pressure  Supply  Line  Pressure  Drop 


FIGURE  5  -  Flow  Capacity  of  Two  Parallel  Regulating  Valves 
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subsonic 


■J  2  2 

m  -  (.00002122)  K2  J|  ?1  ”  ?P 

m  ■  .167  x  2  _  p  2  (2) 

1  P 

where  .00002122  converts  scfh  to  lb/ sec 
Kx  =  834 
K2  -  963 
Cv  *  466 
G  -  1 
T  -  520°R 

x  “  stroke,  2.5  =  full  stroke 

Pj  -  inlet  pressure  =  680  at  design  conditions 

P  =  outlet  pressure  =  400  at  design  conditions 


The  flow  through  the  perforated  plate  may  be  found  with  that  choked  orifice  equation 

AP 


m  -  .287  Pp  (3) 

where  C  *  a  constant,  taken  as  .6 
K  «  .53  for  air 

A  =  .41  (.785)  (8)  *20.6  for  41%  open  holes 
Pp  *  400  at  design  conditions 


(c)  Test  section  flow  -  The  test  section  mass  flow  may  be  found  with  the 
use  of  Chart  1,  NACA  1135  (Ames  Research  Staff,  1953)  or  it  may  be  found,  as  a 
function  of  Mach  number  and  stagnation  conditions,  by  comibining  certain  of  the 
equations  in  the  report  to  give: 
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.  „  A  pstd  Tstd  rmr/.  .  M  .-3 

m  =  po  g  i44  fed  49M(1  +  r> 

M2  1  (5) 

m  -  3.271  PCM  (1  +  j-) 

where  PQ  ■  stagnation  pressure,  psl 
A  =  area,  square  Inches 
To  *  stagnation  temp,  520°R 
M  *  Mach  number 
g  *  gravitational  constant 
pstd  *  .002378  slug/ft3 
Pstd  =  14.7  psia 
Tstd  =  520°R 

At  P0  =  60  psia  and  M  =  1  the  mass  flow  is  113  Ib/sec-  A  nominal  figure  of  115 
lb/sec  is  used. 


The  flow  through  the  subsonic  choke  valve  is  (Jury,  1974): 


m  -  .00002122 


std  „  -  .  f 3417  rsf 

h  c* Pi  sl"  [yr  hj 


since  sin|^3-^-^  J =  1  for  choked  flow 

and  Cg  ■  f  (0)  where  6=valve  position  (Jury,  1974).  Since 
standard  temperature  was  assumed,  then  m  ■  .00002122  P^  f(0). 

Since  the  expansion  of  air  in  a  well  made  nozzle  is  a  highly  efficient 
process,  then  P0=P^  and  T0=T-j-.  Setting  Equations  4  and  5  equal  to  each  other  yidlds 


M2  -3 

3.271  M(1  +y-) 


.00002122  f (0) 
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Equation  6  shows  that  the  test  section  Mach  number  is  a  function  of 
the  valve  position  only. 

(d)  Plenum  flow  -  Treating  the  plenum  flow  passages  as  valves,  the  flow 
coefficient,  Cv,  was  found  to  be  a  function  of  the  seal  stroke,  S,  and  the  stroke 
length,  L,  times  a  constant  experimentally  determined  to  be  26.  There  are  sixteen 
spaces  between  rods,  each  with  a  length  of  about  sixty-four  inches.  The  stroke 

for  maximum  flow  is  approximately  .1  inch.  The  flow,  using  a  subsonic  orifice 
equation  is: 

/  P  2  .  p  2 

m  *  .00002122  K2  Cvl  —  ^ — —  ,  where  Cv  =  26  SL 

For  maximum  total  flow,  S  ■  .1,  L  ■*  16  spaces  at  64  inches.  Using 
G  -  1,  T  -  520°R  and  K2  -  963  the  total  flow  is 

m  -  2.386  /72  I~2 

1  "  2 

^  ^1  taken  as  the  free  stream  static  pressure  of  33.5  psia  (at  60 
psia  PD  and  .95  Mach  number),  and  P2  is  the  pressure  downstream  of  the  subsonic 
choke  valve,  18.2  psia,  then  the  maximum  flow  disregarding  line  losses  is  67  lb/sec. 
However,  since  the  differential  pressure  across  the  rods  at  this  flow  is  much 
greater  than  can  be  tolerated  and  the  flow  requirements  are  only  about  11  lb /sec, 
it  follows  that  the  flow  should  be  throttled.  A  butterfly  control  valve  will  be 
used  for  this  purpose.  This  is  the  plenum  flow  regulating  valve  shown  in  Figure  1. 
The  differential  pressure  required  for  full  flow  is  about  .5  psi  and  all  seals 
opened  fully.  This  is  based  on  a  constant  of  33.5  psia  and  MT0T  of  twelve  per¬ 
cent  of  115  lb/sec  main  line  flow.  The  plenum  flow  desired  is  not  uniform  down 
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the  length  of  the  tunnel.  One  of  the  advantages  of  this  tunnel  is  that  variable 
geometry  can  be  combined  with  a  form  of  variable  porosity.  Therefore  the  seals 

will  not  be  uniformly  opened  and  the  anticipated  differential  pressure  is  1  to  5 
psl. 

3*  — fety  Relief  Flow  Rates.  The  safety  relief  flow  rates  are  based  on 

static  line  conditions  at  maximum  pressure  and  assumed  sudden  malfunction  of 
valves  or  controls  or  system  blockage.  These  flow  rates  are  very  short  tran¬ 
sients  which  decay  rapidly  as  flow  ip  established  and  line  pressures  drop .  They 
represent  worse  conditions  than  can  actually  be  realized. 

(a)  Regulated  pressure  lines  -  The  nominal  line  pressure  is  700  psi. 

The  lines  are  four  inch  schedule  forty  and  eight  inch  schedule  eighty  steel  pipe. 
The  safety  factor  is  four  .at  a  maximum  operating  pressure  of  1157  psi  for  Grade  A 
steel  pipe  (T.O.  00-25-223).  The  lines  Are  projected  by  a  1000  psi  ruptute  disk, 
Figure  1.  The  disk  flow  is  vented  to  atmosphere  above  the  building.  A  steel 
plate  above  the  vent  line  prevents  rain  Water  from  entering  and  prevents  debris 
from  blowing  into  the  air.  An  eight  inch  vent  line  is  used  because  a  four  inch 
line  would  reduce  the  flow  rate  by  either  choking  in  the  vent  or  causing  subsonic 
flow  across  the  ruptured  disk. 

Assume  that  the  PQ  control  Valve  i^  closed  and  that  both  regulators 
fail  to  oben  exposing  the  system  to  the  static  Supply  pressure  of  3000  psi. 

Assume  that  the  disk  ruptures,  leaving  ail  equivalent  of  a  three  inch  round  hole 
with  a  coefficient  of  .8.  Simultaneous  solution  of  the  regulator  and  disk  flow 
equations  shows  that  the  regulator  flow  is  subsonic.1,  the  regulator  and  disk  flow 

is  240  lb/sec,  the  line  pressure  in  the  700  psi  line  is  1830  psi  and  the  safety 
factor  is  2.5. 


(b)  Settling  chamber  and  test  section  -  The  settling  chamber  and  the 
original  test  section  components  were  designed  by  NACA,  Langley  Field,  Virginia 
in  the  1950* s  for  a  pressure  of  45  psig,  and  were  tested  to  a  pressure  of  90 
psig.  New  components  have  been  designed  for  the  same  pressure,  with  safety 
factors  of  four  or  more. 

The  maximum  flow  into  the  settling  chamber  occurs  when  the  regulated 
line  pressure  is  approaching  1000  psi,  the  four  inch  disk  has  not  ruptured  and 
the  PQ  control  valve  fails  wide  open.  Simultaneous  solution  of  the  control  valve 
and  perforated  plate  flow  equations  shows  that  the  pressure  upstream  of  the  per¬ 
forated  plate  is  810  psi  and  the  flow  is  225  lb/sec.  The  worst  case  is  when  the 

subsonic  choke  valve  is  also  closed  and  all  the  flow  must  be  exhausted  through  re¬ 
s'  : 

lief  devices.  Two  such  devices  are  in  use.  One  is  a  relief  valve  set  at  65  psia. 
The  valve  has  an  eight  inch  inlet,  twenty-six  square  inch  flow  area,  ten  inch 
outlet  and  expands  into  a  twelve  inch  pipe  which  routes  the  flow  through  the 
tunnel  exhaust  and  silencer.  The  second  is  a  twelve  inch,  70  psia  rupture  disk. 

An  iterative  solution  shows  that  the  pressure  rises  to  120  psig  for  a  safety 
factor  of  1.5.  Flow  through  the  ruptured  disk  is  subsonic  at  175  lb/sec.  Flow 
through  the  relief  valve  is  sonic  at  50  Ib/sec.  The  safety  factor  is  inadequate, 
however,  it  was  calculated  for  a  line  pressure  of  300  psi  above  the  set  point  and 
for  simultaneous  failure  of  two  valves.  Interlocks  will  be  used  to  increase  the 
safety  of  the  facility. 

(c)  Interlock  system  -  The  interlock  system  is  designed  to  shut  off  the 
high  pressure  air  supply,  and  hence  shut  the  facility  down,  in  the  event  of: 
overpressure  of  the  regulated  pressure  lines  or  settling  chamber,  operation  of 
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relief  devices,  hydraulic  system  failure,  or  the  choke  valve  closed  beyond  normal 
operation.  The  devices  chosen  to  shut  the  facility  down  in  an  emergency  are  the 
dome  loaded  regulating  valves  and  the  isolation  valve  because  they  can  be  made  to 
close  in  the  absence  of  electrical  or  hydraulic  power.  The  isolation  valve  is 
operated  via  a  normally  closed  pilot  valve  and  air  pressure.  The  domes  of  the 
regulators  are  pressurized  with  an  electrically  operated  dome  loading  regulator 
which  bleeds  off  the  dome  pressure  with  no  electric  power  applied. 

In  the  design  of  the  system,  consideration  was  given  to  the  fact  that 
in  protecting  against  one  hazard,  another  hazard  can  be  created.  An  example 
is  possible  damage  to  the  model  if  the  model  is  not  at  zero  degrees  pitch  and  if 
the  tunnel  is  running  supersonically  and  goes  out  of  flow.  This  occurs  when  the 
normal  shock  passes  through  the  test  section.  Because  of  this  hazard  to  the 
model,  the  interlocks  for  most  events  incorporate  a  short  time  delay  to  allow 
the  model  to  pitch  to  zero  and  also  start  the  pitch  motion.  Closing  of  the  choke 
valve  however  initiates  immediate  tunnel  shutdown  because,  as  seen  in  Figure  8, 
the  tunnel  will  run  subsonically  with  the  valve  less  than  fifty-three  degrees  open. 

B.  HYDRAULIC  VALVE  POSITIONING  SYSTEMS 

Hydraulic  power  was  selected  for  valve  positioning  because  of  the  need  for 
high  performance  in  a  short  running  time  blowdown  tunnel.  Two  valves,  the  stagna¬ 
tion  pressure  control  valve  and  the  subsonic  choke  control  valve,  will  be  hydrau¬ 
lically  operated.  A  third  valve,  the  plenum  chamber  suction  control  valve,  may  be 
driven  in  the  future  although  initially  the  valve  will  be  manually  preset. 

The  stagnation  pressure,  PQ,  control  valve  will  be  initially  set  up  for  closed 
loop  valve  position  control.  This  is  the  checkout  and  pre-regulator  setting 
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purposes.  The  subsonic  choke  valve  will  be  closed  loop  position  controlled  for 
purposes  of  open  loop  control  of  subsonic  Mach  number  based  on  prior  tunnel 
calibration  data. 


A  block  diagram  for  the  servo  valve  and  load  is  shown  in  Figure  6.  A  Moog 
Series  73  valve,  rated  at  5  gpm  is  used  for  each  system.  The  servo  valve  second 
order  model,  shown  in  Figure  6,  is  accurate  to  about  300  radians/second  with  a 
natural  frequency,  wgv,  880  radians  per  second  and  a  damping  ratio,  tOgv* 

0.9  (Moog  Catalog).  The  flow  gain  is  Kj.,  with  units  of  cubic  inches  per  second 
(cis)  per  ma,  and  is  discussed  as  a  part  of  load  dynamics. 


The  viscous  damping  factor,  B,  usually  low,  is  assumed  to  be  zero.  Applying 

block  diagram  reduction  techniques  to  the  load  portion  of  Figure  6  gives  a  trans- 

i 

fer  function  of 

1/A _ _ 


fw 


+  !£*  +  +  ^ 

A2  A2  A2  A 2 


where  x  “  Piston  travel,  inches  (load  travel**piston  travel  or  is  referred 
to  piston  travel) 

Q  -  Load  flow,  cis 
A  =  Net  piston  area,  square  inches 

K3«  Compressibility  coefficient,  V/B,  in  cubic  inches/psi 
(Johnson  and  Schmid,  1956) 

V  m  Volume  under  compression, 

B  =  Bulk  modulus,  170,000  psi  for  mineral  based  oil 
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2 

M  *  Load  mass,  lb/inch/sec 

K  =  Load  spring,  lb/inch 
s 

1^“  Pressure  gain,  cis/psi  (often  referred  to  as  C^,  with 
C2  being  a  laminar  leakage  coefficient) 

The  arrangement  of  the  transfer  function  provides  a  single  equation  for  use 
with  either  a  pure  mass  load  or  for  a  load  with  a  parallel  spring.  The  low  fre¬ 
quency  behavior,  u)  approaching  zero,  is  that  of  an  integrator  or  a  gain  of 
2 

For  most  systems  A  ,  is  much  greater  than  K  K~  and  at  high  frequency  the  constant 

S  j 

term  in  the  denominator  is  negligible  compared  to  the  other  terms.  Then  at  high 
frequencies  the  denominator  is  the  time  constant  form  of  a  second  order  system 
with  a  free  integration  whose  natural  frequency  and  damping  ratio  are  given  by: 


The  constants  A,  M  and  Kg  are  different  for  each  system  and  are  discussed 
separately.  The  constants  and  are  taken  froiff  a  graph  of  the  servo  valve 
pressure,  flow  and  current  characteristics.  Figure  7. 


FIGURE  7  -  Characteristics  of  5  gpm  Servovalve  with 
3000  psi  Supply  Pressure 
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Both  systems  are  lightly  loaded,  operating  in  the  linear  region  at  the  left  of 

the  graph.  The  flow  gain,  1^,  is  nearly  constant  in  this  region  at  2.13  cis/ma. 

The  pressure  gain,  K2>  is  not  constant,  varying  from  nearly  zero  at  the  origin  to 
-3 

5.85  x  10  cxs/psi  at  the  top  of  the  graph,  marked  "Operating  Region".  Hansen 
(1961)  gives  a  typical  order  of  magnitude  value  of  10~3  cis/psi.  Therefore  a 
median  value  of  3  x  10  “3  is  used. 


1*  -P-ressure  Control  Valve.  The  load  mass,  a  linear  motion  air  valve  plug, 
is  .16  lbs/in/sec2.  The  actuator  is  a  2%  inch  bore,  1  inch  double  ended  rod 
cylinder  with  an  area  of  4.12  square  inches.  The  stroke  is  2.5  inches.  With 
the  piston  centered,  the  volume  under  compression  on  one  side  is  5.15  cubic 
inches.  The  constant  resisting  force,  whence  Ks»0,  is  caused  by  differential 
air  pressure  acting  over  a  small  unbalanced  area  on  the  virtually  balanced  double 
acting  air  control  valve  plug  and  is  800  pounds.  The  load  pressure  drop  is  800/ 
4.12  =  194  psi.  The  flow  at  194  psi  and  15  ma,  Figure  7,  is  31  cis  for  a 

maximum  velocity  of  31/4.12  *  7.52  in/sec.  The  pressure  drop  across  the  valve  is 
3000  -  194  =  2806  psi. 

The  compressibility  coefficient,  is  3  x  10"5.  The  natural  frequency, 
Wl,  is  1880  radians  per  second.  Multiplying  the  servo  valve  and  load  transfer 
functions  together  would  produce  a  characteristic  equation  whose  roots  would  show 
a  frequency  of  880  radians/second.  In  other  words,  the  servo  valve  cannot  control 
a  mass  system  at  that  high  a  frequency.  The  performance  must  be  downgraded  so  that 
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attenuation  occurs  before  u)sv.  Then  the  high  frequency  load  characteristics  are 
reduced  to  almost  nothing  and  the  load  equation  can  be  simplified  to: 


Crossover  should  occur  at  about  1/3  (x)gv  (Broome,  1963)  or  300  radians/ 
second*  With  no  poles  or  zeros  to  the  left  of  0)c,  then  toc  «  Kv,  the  velocity 
coefficient*  The  gain  will  be  adjusted  to  achieve  full  velocity.  With  a  feed¬ 
back  coefficient,  Kfb,  of  six  volts/inch,  the  amplifier  gain  is: 


(Ky)  (imax)  (300)  (15) 
K  "  (Kfb) (Vel)  "  (6) (7.5) 


100  ma/volt 


This  is  a  well  behaved  second  order  system  with  a  free  integration.  Itfs 
relative  stability  can  be  readily  shown  to  be  adequate.  Slowing  the  system  down 
to  match  the  servo  valve  is  only  a  degration  of  performance  in  the  sense  of  fol¬ 
lowing  an  oscillatory  input.  It  will  be  shown  in  Section  iII-C  that  the  cross¬ 
over  frequency  will  be  further  reduced  but  that  a  high  velocity  is  desirable  and 
is  retained. 


2.  Subsonic  Choke  Valve.  The  load  mass  is  a  rotary  disk  in  a  butterfly  style 
valve  which  is  driven  by  a  crank,  and  linear  hydraulic  cylinder.  Operating 
characteristics  of  the  valve  are  shown  in  Figure  8,  Curves  A  through  G. 

Curve  A  shows  the  rated  flow  through  the  valve.  Curve  B  shows  the  Mach 
number  at  various  disk  angles.  Past  the  angle  for  Mach  1  the  flow  cannot  be 
choked  in  the  valve  and  the  valve  flow  will  be  subsonic.  Tunnel  Mach  number  will 
then  depend  on  other  conditions. 
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Curve  C  shows  the  typical  unit  torque  characteristics  (Beard,  1969)  of 
butterfly  valves.  The  high  torque  peak  at  eighty  degrees  will  not  be  experi¬ 
enced  by  the  valve  because  the  differential  pressure  drops  as  shown  on  Curve  D 
when  the  valve  unchokes. 

Curve  E  shows  the  stroke  versus  valve  disk  angle.  The  feedback  potentiom¬ 
eter  is  rotary,  which  introduces  a  slight  nonlinearity  in  the  analysis.  The 
coefficient  is  selected  at  the  design  point. 

Curve  F  shows  the  unit  force  per  unit  torque  which  when  multiplied  by  the 
values  of  Curves  C  and  D  provides  the  operating  force  curve  shown  on  Curve  G. 

The  pi6ton  stroke,  crank  length  and  initial  angle  were  chosen  for  minimum  F/T 
and  maximum  linearity  of  F/T  in  the  Mach  number  range  of  .3  to  .9.  The  operating 
forces  to  move  the  disk  shown  on  Curve  G  provide  a  variable  mechanical  spring 
constant.  Since  the  forces  are  caused  by  airflow,  and  the  valve  will  be  operated 
with  and  without  flow,  the  spring  constant  varies  from  zero  to  a  maximum  at  the 
design  point. 

Performance  required  of  the  control  system,  other  than  stability,  is  mainly 
accuracy.  As  Mach  1  is  approached,  the  curve  of  Mach  number  versus  disk  angle 
gets  very  steep  as  shown  on  Curve  B.  The  valve  is  normally  pre-set  before  a  run 
with  possibly  minor  adjustments  made  during  a  run.  Hence  operating  speed  is  not 
a  problem. 

(a)  Design  analysis  -  The  valve  disk  is  approximated  by  a  steel  six  inch 

diameter  hub,  sixteen  inch  long  plus  a  one  inch  thick,  fifteen  inch  diameter  disk, 

2 

having  a  combined  polar  moment  of  inertia,  J,  of  3.20  inches  lb  sec  .  The  linear 
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motion  diagram  of  Figure  6,  and  the  resulting  load  transfer  function,  apply  to 
this  system  if  the  rotating  components  are  referred  to  the  linear.  The  rotat¬ 
ing  mass  can  be  referred  to  the  linear  piston  operation  by: 

JB  -  Mrx 

where  r  *  radius  «  from  Figure  8,  Curve  F 
F/T  ■  .14  at  design  point  and 
x  -  .1259 

since  x  “  .1256  at  design  point. 

Then  M  -  3.58  lb/in/sec2 

The  hydraulic  cylinder  has  a  2h  inch  bore,  a  1  inch  rod  and  is 
pivoted  at  one  end.  The  average  area,  A,  is  4.51  square  inches.  Using  the 
average  area  and  the  design  load  resisting  force,  2800  lb,  the  load  pressure 
drop  is  621  psi.  Then  at  maximum  current,  from  Figure  7  the  flow,  Q,  is  28.2 
cis  and  the  velocity,  Q/A,  is  6.25  in/sec.  This  is  an  approximation  with  an 
unbalanced  piston  area.  The  velocity  can  be  calculated  from  a  balance  of  forces. 
The  valve  is  rated  by  Moog  at  5  gpm  at  1000  psid  valve  drop.  The  drop  through 
each  side  of  the  valve  is  then 


AP  -  500  J  —'\\7  Cl^~ 

*  AV 

where  A  is  area,  ih2  and  V  is  velocity,  in/sec.  Using  3000  psi  sup 
Ply»  Ag  for  the  blind  end  area  and  AR  for  the  rod  end  area,  a  balance  of  forces 
gives: 
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The  velocities  found  by  the  two  different  methods  differ  very  little 
because  the  ratio  of  rod  area  to  blind  area  is  .84  or  nearly  unity.  An  average 
piston  area  will  be  used  and  the  areas  assumed  equal  to  avoid  complicating  the 

analysis. 

At  mid  stroke  the  volume  of  oil  under  compression  on  one  side  is  22.51 
cubic  inches,  neglecting  the  rigid  connecting  tubing.  The  compressibility  coef¬ 
ficient,  K3,  is  then  1.33  x  lO-4.  The  natural  frequency,  wL,  is  207  radians  per 
second.  Using  the  same  rule  of  thumb  by  Broome,  crossover  occurring  at  about  1/3 
WL,  wc  should  be  about  70  radians  per  second.  With  an  attenuation  of  18  DB  per 
octave  after  207  radians  per  second,  the  attenuation  is  so  great  at  880  radians 
per  second  that  the  servo  valve  characteristics  can  be  ignored.  The  load  dynamics 
dominate  the  performance. 

The  damping  ratio,  «*,,  is  .054.  Using  a  small  tube  to  bypass  fluid 
from  the  piston  inlet  to  exhaust  can  improve  the  damping  ratio.  A  flow  of  2  gpm, 
or  7.69  cis  at  1000  psid  can  be  tolerated.  Since  the  load  pressure  drop  is  621 
psid,  the  maximum  quiescent  flow,  under  load,  is  1.24  gpm.  A  10  gpm  pump  is 
available  for  supplying  leakage  flow  plus  a  small  accumulator  for  supplying  peak 
demand  of  17  gpm  for  less  than  1.5  seconds. 

At  7.69  cis  per  1000  psid,  the  crossport  leakage  coefficient  is  7.69  x 
10-3  Cis/psi.  This  improves  K2  from  3  x  10_3  to  10.69  x  lO-3.  The  damping  ratio 
increases  from  0.54  to  .192. 

The  feedback  coefficient,  Kfb,  is  1.33  volts/inch  at  the  operating 
point,  the  gain  to  operate  at  a  Ky  -  u>c  -  70,  with  no  poles  or  zeros  to  the  left 
of  wc,  and  at  the  calculated  velocity  of  6.25  inches/second  is: 


K 


126  ma/volt 


(Kv>(1mrt 

®fb)<T«x) 


The  steady  state  gain  with  a  spring  constant  of  1800  lb/in  is: 


“Ab 


K2Ks 


-  137 


To  provide  separate  adjustment  of  high  and  low  frequency  gain,  a  lag 
compensator  is  added  and  the  gain  is  increased  by  75Z.  The  amplifier  gain  becomes 
220  ma/volt,  the  loop  gain  approximately  140  and  the  steady  state  gain  with  1800 
lb/in  spring  constant  of  146.  The  compensator  provides  a  gain  margin  at  (o^of  2. 
Additional  gain  margin  is  provided  by  the  neglected  and  somewhat  nebulous  load 
damping.  The  transfer  function  for  the  lag  compensator  is: 


v0  R2CS*1 
vj  S  (Ri  +  R2)CS*1 

Where  C  =  5uf 

R]  =  236Kft 
R2  =  .50101 


An  open  loop  frequency  response  plot  of  the  system  is  shown  in  Figure 
9  with  spring  constants  of  0,  500  and  2000  lb/ln.  Also  shown,  labeled  "before". 
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is  the  high  frequency  portion  of  the  system  prior  to  the  addition  of  the  lag  filter, 
gain  increase  and  crossport  leakage.  The  transfer  function  plotted  is: 

vo  (T.S+1) 

^(s)  ’  K  C^S+IT  (2.13)^--=- 

- - - — - — — - - (i.33) 

(2 . 34x10~5S3+.  176K2S2+  (1+6 . 54xl(f  \)  ST .  0491^ 


where  the  values  used  are: 


K 

T, 


K 


INITIAL 


126 

0 

0 

3x10 
0 


-3 


FINAL  SYSTEM 


220 
125 
1.43 
10.79x10 


“4 


0,  500,  2Q00 


The  crossport  leakage  decreases  the  maximum  velocity  under  load  from 
6.25  in/sec  to  5.19  in/sec.  The  leakage  also  dissipates  .73  horsepower  under  load. 
The  heat  dissipated  results  in  a  gradual  oil  temperature  buildup  which  is  not  a 
problem  fpr  a  blowdown  tunnel. 

Accuracy  is  difficult  to  address  frcfcn  theory.  First  it  must  be 
assumed  that  the  two  pivots  have  minimal  clearances  and  that  the  typically  bother¬ 
some  keyway  has  no  clearance,  possibly  a  light  drive  fit  on  the  sides.  The  possi¬ 
bility  of  loss  motion  between  the  positioning  cylinder  and  the  measuring  point  can 

1 

introduce  other  control  problems  besides  accuracy.  With  no  lost  motion,  there  are 
several  ways  to  attempt  to  assess  accuracy. 
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AMPLITUDE  RATIO 


FIGURE  9  -  Butterfly  Choke  Valve  Frequency  Response 


Thfs  steady  state  characteristics  are  that  of  a  proportional  control  sys¬ 
tem  when  tunnel  air  is  flowing.  Proportional  control  has  an  error  of  the  step 
command  divided  by  one  plus  the  loop  gain.  If  the  step  is  6.3  inches  then  the 
error  is  .033  inches.  However  it  is  probable  that  only  the  last  half  inch  is 
slow  enough  to  be  considered  a  first  order  response.  In  any  event,  the  .033 
inches  error  produces  7.7  ma,  half  scale,  from  the  amplifier,  which  is  certainly 

going  to  move  the  load. 

If  the  disk  is  assumed  to  be  exactly  positioned  and  then  the  load  is 
applied,  the  piston  will  have  to  move  .02  inch  to  correct  for  oil  compression. 
The  proportional  droop  is  then  .0001  inch. 

If  it  requires  a  3Z  of  maximum  signal  to  move  the  servo  valve  (an  assump¬ 
tion  of  Broome),  then  the  error  which  can  be  sustained  is  .002  inch  with  the 
amplifier  and  feedback  gains  is  use.  This  represents  one  minute  of  angle  on 
the  disk.  It  is  difficult  to  read  a  precision  inclinometer  closer  than  two 
minutes  of  angle.  Therefore  the  accuracy  is  probably  as  good  as  can  be  rea¬ 
sonably  put  to  use. 

C.  PRESSURE  CONTROL  SYSTEM 

1,  Description  and  Coefficients 

(a)  Description  of  system  -  This  is  a  dual  volume  system  which  is 
modulated  by  the  PQ  control  valve.  Position  control  of  the  PQ  valve  for  check¬ 
out  purposes  was  discussed  under  Section  II-B.  A  schematic  diagram  and  a 
block  diagram  of  the  system  is  shown  in  Figure  10  and  11,  respectively.  The 
constants  used  in  the  figures  are  discussed  later  and  given  in  Table  I. 
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It  is  assumed  that  the  flow  rate  is  sufficiently  high  to  cause 
choked  flow  at  the  perforated  plate  and  at  either  the  nozzle  throat  or  the 
subsonic  choke  valve.  The  choice  of  the  latter  choke  point  depends  on  Mach 
number  and  not  pressure.  This  assumption  is  valid  when  stagnation  pressure,  PQ, 
is  above  30  psia  since  the  exhaust  is  to  atmosphere.  The  pressure  upstream  of 
the  perforated  plate,  Pp,  must  be  greater  than  60  psia  since  the  exhaust  is  to 
PD.  The  justification  for  this  assumption  is  twofold.  Firstly,  the  choked 
flow  at  Pp  of  60  psia  is  only  17  lb/sec  or  15%  of  maximum.  This  is  a  starting 
condition  of  little  interest.  Secondly,  a  negative  feedback  coefficient  which 
is  zero  based  on  the  assumption,  is  non-zero  when  the  assumption  is  invalid. 

Since  negative  feedback  has  a  stabilizing  effect,  the  system  stability  increases 
when  the  assumption  is  invalid. 

A  perturbation  analysis  is  required  due  to  nonlinearities.  At  Mach 
1  or  above  the  system  is  inherently  more  stable  due  to  the  higher  mass  flows 
which  produce  a  greater  negative  feedback.  However  at  this  time,  with  60  psia 
controlled  pressure,  the  control  valve  requires  a  higher  velocity  to  stay  within 
the  linear  model  because  the  valve  is  not  choked  and  there  is  a  downstream  pres¬ 
sure  effect.  At  Mach  .3,  the  lowest  Mach  number  of  interest,  the  reverse  is  true. 
Inherent  stability  is  less  but  the  valve  velocity  requirements  are  less  and  the 
large  power  application  which  can  be  demanded  by  phase  lead  type  of  compensation 
can  be  realized. 

(b)  System  coefficients  and  characteristics  -  The  constants  and  coef¬ 
ficients  which  apply  to  Figure  11  are  tabulated  in  Table  I.  Some  of  the  coeffici¬ 
ents  are  partial  derivatives  of  equations  numbered  in  Section  III-B,  with  respect 
to  the  variable  listed  in  the  table. 
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TABLE  I 


CONSTANTS  AND  COEFFICIENTS  FOR  PRESSURE  CONTROL 


APPLICABLE 
FOR  MACH 

QUANTITY 

VALUE 

NUMBER  OF 

SOURCE  AND  COMMENTS 

ha 

695  psi 

.3 

Constant 

p1b 

669  psi 

Constant 

R 

12  x  53-3 

.3,  1 

Constant 

TP 

520°R 

•3,  1 

Constant 

To 

520°R 

•3,  1 

Constant 

VP 

3317  cu  in 

•  3,  1 

Constant 

VOA 

215,800  cu  in 

.3 

Constant 

V0B 

227,260  cu  in 

1 

Constant 

k2A 

96.3  lbs/sec/?n 

•  3 

Equation  1 ,  Var  X, 

Son i c  Valve  FI ow 

k2B 

86.2  lbs/sec/ in 

1 

Equation  2,  Var  X, 

Subsonic  Valve  Flow 

k3A 

0 

.3 

Equation  1,  Var  Pp, 

Sonic  Valve  Flow 

k3B 

-.151  ibs/sec/psi 

1 

Equation  2,  Var  Pp, 

Subsonic  Valve  Flow 

*A 

.287  Ibs/sec/psi 

•3,  1 

Equation  3,  Var  Pp, 

Sonic  Perforated  Plate  Flow 

k5 

0 

.3,  1 

Equation  3,  Var  P0, 

Sonic  Perforated  Plate  Flow 

K6A 

.933  Ibs/sec/psi 

•3 

Equation  A,  Var  PQ, 

Sonic  Choke  Valve  Flow  ■ 

Subsonic  Nozzle  Flow 

k6b 

1.927  Ibs/sec/psi 

1 

Equation  A,  Var  P0, 

Sonic  Nozzle  Flow 

The  system  of  Figure  11  contains  two  first  order  systems  in  cascade.  Apply¬ 
ing  block  diagram  algebra  to  the  first  system,  for  variable  x  to  variable  Pp,  and 
using  the  A  subscripts,  results  in  a  transfer  function  of  the  form: 
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K, 


2A 


K4  -  K3A 


K 


VP  s 
RTp  (k4-k3A) 


+  1 


£+1 

0) 


where  K  ■  steady  state  gain 
a)  -  break  frequency 

Using  this  form  of  the  transfer  function,  and  the  values  given  in 
Table  I,  the  characteristics  have  been  computed  and  are  given  in  Table  II. 


TABLE  II 

PRESSURE  SYSTEM  CHARACTERISTICS 


NAME 

MACH  NUMBER 

GAIN 

BREAK  FREQUENCY 

Flow  into  Vp 

.3 

335.5  psi/in 

28.78 

Flow  into  Vc 

•  3 

.308  psi/psi 

1.438 

Flow  into  Vp 

1 

196.8  psi/in 

43.92 

Flow  into  VQ 

1 

. 149  psi/psi 

2.820 

(c)  Valve  positioning  coefficients  -  It  was  shown  in  Section  II-B,  that 
the  load  dynamics  of  driving  the  valve  plug  mass  may  be  simplified  to  l/SA,  where 
A  is  4.12  square  inches.  Similarly,  based  on  the  break  frequencies  given  above, 
the  output  signal  will  be  greatly  attenuated  before  the  servo  valve  natural  fre¬ 
quency,  b)gv,  is  reached.  Therefore  the  servo  valve  may  be  approximated  by  a  pure 
gain  of  2.13  cis/ma. 
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(d)  Performance  requirements  -  A  high  degree  of  accuracy  is  required  so 
that  tunnel  conditions  are  repeatable.  This  can  help  in  troubleshooting  data 
problems  using  raw  data.  Reduced  data  is  typically  non-dimensionalized. 

The  maximum  steady  state  airflow  of  the  tunnel  is  more  than  fifty 
times  the  rate  of  which  the  storage  cylinders  can  be  recharged.  This  means  that 
once  an  appreciable  mass  flow  has  started  the  tunnel  should  be  on  condition  in  a 
minimum  time.  In  addition,  valve  opening  and  closing  should  be  smooth  and  should 
require  a  couple  seconds  to  avoid  water  hammer  in  the  supply  lines.  This  require¬ 
ment  is  translated  to  mean  three  to  four  seconds  from  start  to  operating  condition 
and  2.5  second  shutdown.  Time  at  operating  condition  may  be  as  little  as  ten 
seconds  if  model  attitude  and  tunnel  walls  are  not  changed,  or  it  may  be  much 
longer . 

Operation  should  be  simple.  It  is  preferred,  but  not  required,  that 
no  controller  adjustments  are  required  in  changing  from  one  Mach  number  to  another. 

2.  Frequency  Domain  Design.  The  system  is  partially  compensated  by  feeding  the 
intermediate  pressure,  Pp,  back  through  a  high  pass  filter  in  a  minor  loop.  This 
is  a  phase  lead  type  of  compensator  which  has  several  advantages  over  phase  lead 
compensation  in  the  forward  path.  It  provides  a  fixed  zero  but  a  pole  which 
increases  in  frequency  as  gain  is  increased.  It  provides  a  means  of  independently 
setting  the  gain  margin,  using  the  error  amplifier  gain,  and  the  phase  margin, 
using  the  minor  loop  feedback  gain.  It  also  provides  a  means  of  dampening  oscil¬ 
lation  in  Pp. 
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Lag  compensation  is  also  used  to  permit  a  gain  increase  and  improve 
the  steady  state  accuracy.  Theoretically  this  should  not  be  required  because  the 
free  integration  provided  by  the  hydraulic  cylinder  provides  zero  steady  state 
error.  In  a  real  system  however,  the  presence  of  unknown  forces  such  as  stiction 
can  modify  the  integral  effect  at  steady  state  to  that  of  a  gain.  The  accuracy 
requirement  and  the  servo-valve  performance  determine  the  gain.  The  strain  gage 
pressure  transducer,  used  for  feedback,  is  accurate  to  about . .25%  at  60  psia. 

The  servo  system  will  not  degrade  the  performance  if  it  has  an  accuracy  of  .15% 
of  60  psia,  or  .0225  volts  at  .25  volts  per  psi.  If  the  servo  valve  requires  3% 
of  the  15  ma  full  scale  to  move,  then  .0225  volts  into  the  error  amplifier  should 
provide  .45  ma  out.  The  required  gain  is  20  ma/volt.  A  gain  of  44  will  be  used 
at  Mach  1. 

The  lag  compensator  is  not  intended  to  provide  an  accurate  ramp  fol¬ 
lowing  capability  because  the  large  capacitor  required  would  make  the  system 
sluggish  in  settling  down  to  a  steady  state  value.  A  ramp  input  will  be  used, 
however  an  input  prefilter  will  be  used  to  round  the  corners  of  the  ramp  and  pro¬ 
vide  smoother  operation  and  decreased  settling  time  at  the  operating  pressure. 

A  condensed  block  diagram  of  the  system  is  shown  in  Figure  12.  Addi¬ 
tional  system  gains  and  pole  locations  for  Mach  .1  operation  are  also  shown.  The 
minor  loop  gains  and  feedback  filter  were  selected  through  the  use  of  a  Nyquist 
diagram  plotting  program.  Values  for  the  Mach  1  case  were  adjusted  until  the 
minimum  values  of  both  gain  margin  and  phase  margin  as  suggested  by  Doebelin  (1962) 
were  achieved.  These  are  2.5  and  30  degrees  respectively.  The  forward  path  gains, 
Ka  and  Kg  on  Figure  12,  were  then  adjusted  until  adequate  response  was  achieved. 
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(See  "Time  Domain  Simulation").  Nyquist  plot  of  the  open  loop  frequency  response 
of  the  system  shown  in  Figure  12,  with  the  attenuator  set  at  1,  is  shown  in 
Figure  13A. 


Attenuator ' 
Lag 

compensator 


Servo  valve 

Hydraulic  cylinder  &  Test  section 
air  valve  flow  airflow 


m 

A 

B 

c 

D 

.i 

346.2 

28.78 

.883 

.501 

.525 

.3 

335.5 

28.78 

* 

O 

00 

1.438 

.619 

1 

196-8 

43-92 

.149 

2.820 

1 

FIGURE  12  -  PRESSURE  CONTROL  SYSTEM-CONDENSED  BLOCK  DIAGRAM 


The  attenuator  of  Figure  12  is  a  potentiometer  which  is  geared  to  the 
choker  valve  shaft.  The  potentiometer  has  special  taps  at  the  approximate  angles 
for  Mach  .1,  .3  and  1.0  operation  of  the  choker.  Shunting  resistors  will  be  used 
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FIGURE  13  -  Pressure  Control  System  Frequency  Response  -  Nyquist  Plot 


across  these  taps  to  provide  the  attenuation  shown  in  Figure  12  for  MN  .1  and  .3 

This  corrects  the  gain  and  phase  margins  for  low  Mach  number  operation  as  shown 

in  Figure  13B.  The  gain  is  then  automatically  adjusted  and  no  operator  adjust- 
ment  is  required. 

An  open  loop  Bode  plot,  corresponding  to  Figure  13B,  is  shown  in 
Figure  14.  The  lag  compensator  pole  to  zero  frequency  ratio  was  chosen  to  pro¬ 
vide  the  minimum  gain.  The  pole  was  selected  at  a  fairly  high  frequency  which 

would  still  give  the  desired  effect  and  would  keep  the  capacitor  and  resistors 
reasonably  sized. 

3*  — —  P0™1^.. Simulation.  The  system  was  simulated  using  an  IBM  1130  Con¬ 

tinuous  Systems  Modeling  Program  which  runs  on  the  9"  x  9"  Facility  Hewlett- 
Packard  mini-computer. 

For  the  simulation,  the  gain  blocks  K2  and  K3  shown  in  Figure  11  were 
replaced  with  a  subroutine  which  simulated  the  actual  Pc  control  valve  operation. 
The  subroutine  considered  upstream  line  losses  and  calculated  flow  for  either  the 
sonic  or  subsonic  case,  depending  on  the  pressure  drop  available.  A  block  diagram 
of  the  subroutine  is  shown  in  Figure  15. 

A  ramp  rate  of  22.5  psi/sec  was  used.  This  is  a  two  second  ramp  from 
ambient  to  maximum  pressure  of  45  psig.  Since  simulation  runs  were  started  at 
15  psig  only  1.33  seconds  was  required  for  the  ramp. 
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FIGURE  14  -  Pressure  Control  System  Frequency  Response  -  Bode  Plot 


Where: 

me  =  Estimated  Flow 
m0  s  Actual  Flow 
P|  =  f(mQ)  due  to  line  loss 
Pp/Pj  =  Pressure  Ratio,  .53=  Sonic 
E  =  Flow  estimation  error 
Pr  s  Reguloted  Pressure  s  700psia 


I 


Control 

Valve 

Perforated 

Plate 


FIGURE  15  -  Block  Diagram  of  Control  Valve  Subroutine 


The  simulation  runs  showed  that  the  state  variables  were  never  saturated 
and  that  the  linear  model  was  a  good  approximation.  The  maximum  valve  velocity 
and  position  of  7.5  in/ sec  and  2.5  inches  respectively  were  never  reached. 

Several  runs  were  made  with  different  input  prefilter  time  constants  at 
Mach  1.  The  shortest  time  to  0.5%  settling,  2.20  seconds,  was  obtained  with  a 
prefilter  time  constant  of  .3  seconds.  The  additional  time  required  for  starting 
at  0  psig  instead  of  15  psig  is  .67  seconds  for  a  total  time  of  less  than  3  seconds. 
Since  saturation  did  not  occur,  the  ramp  rate  can  be  increased  to  30  psi/sec  or 
more.  The  total  time  to  stabilized  conditions  is  less  than  2.5  seconds  and  the 
air  lost  in  the  process  is  less  than  70  pounds  mass. 

The  velocity  constant,  Ky,  recorded  was  17.7  which  is  a  little  less  than 
the  theoretical  value  for  the  linear  system.  Four  of  the  simulation  runs  for 
Mach  1  are  shown  in  Figure  16.  Values  recorded  for  the  runs,  including  variables 
subject  to  saturating,  are  shown  in  Table  III. 


TABLE  III 


DATA  RECORDED  IN  Pc  SIMULATION  RUNS 


ITEM 

VALUE 

RUN 

A 

B 

c 

D 

Pre-filter  Time  Constant 

.01 

.25 

.3 

.4 

Maximum  Control  Valve  Speed 

-2.27 

+1.09 

-2.81 

+.90 

Maximum  Valve  Position 

1.61 

1.44 

1.40 

1.32 

Maximum  Flow 

129.8 

122.6 

120.9 

117.0 

|| 

Maximum  P0 

61.79 

60.33 

60.29 

60.23 

0. 5%  Setti ing  Time 

3-55 

3.10 

2.20 

2.55 

TIME,  SEC 


FIGURE  16  -  Pc  Simulation  Runs-Rosponsc  to  a 
Rarap  to  Operating  Pressure 


VI'.  CONCLUSIONS 


The  existence  of  a  design  source  document  on  the  wind  tunnel  should 
prove  valuable  in  the  future  for  the  familiarization  of  personnel  for 
maintenance,  operation  and  any  future  modifications  of  the  tunnel.  De¬ 
tailed  design  information  is  given  which  includes,  for  instance,  the 
trends  in  airflow  control  systems  stability  associated  with  changing 
tunnel  Mach  number.  The  piping  and  control  systems  are  shown  to  be 
adequate  from  the  standpoints  of  safety,  control  accuracy  and  facility 
running  time. 
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